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Abstract 

In this article, one defines two models of adaptive systems: the meta- dynamical 
adaptive system using the notion of Kalman dynamical systems and the adaptive 
differential equations using the notion of variable dimension spaces. This concept of 
variable dimension spaces relates the notion of spaces to the notion of dimensions. 
First, a computational model of the Douady's Rabbit fractal is obtained by using 
the meta-dynamical adaptive system concept. Then, we focus on a defense-attack 
biological model described by our two formalisms. 
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1 Introduction 



In the two last decades, there has been much interest in the study and for- 
malization of complex adaptive systems (see [21] and [S]). Many different ap- 
proaches have been proposed: artificial chemistries, evolutionary formalism, 
cellular automata. However, from a more theoretical viewpoint, few mathe- 
matical formalisms exist for adaptive systems. Though we may cite the chap- 
ter "Categorical System Theory" proposed by A.H. Louie in [18] who discuss 
the relationship between natural and formal systems, most of the studies are 
mainly simulation-based (see [8] and [9]). With the recent evolution of physics 
and biology, a general mathematical formalism for adaptive systems would be 
very useful. Actually, behind their apparent heterogeneity, most of the adap- 
tive systems share one important feature: they are dynamical objects whose 
structures are sometimes modified by a top level automation-like rule. On the 
basis of this observation, we have built a two-level formalism that helps us to 
design an algorithm for fractals and a biological model of co-evolution in a 
bacterium-phage system. 

The main goal of this article is to present a concept of adaptive systems, gen- 
eral enough to be used in different fields (mathematics, physics or biology) and 
to apply this framework to a fractal implementation and a biological system. 
Besides, a new concept of spaces is developed. Since the second part of the 
19*^ century, various definitions of dimensions, in particular those developed 
by Cantor and Peano, have appeared (see [15 )• These new definitions lead to 
the concept of a fractal model and with it, a global view of spaces (see [13] )• In 
order to study adaptive systems, we will use a sort of "adaptive" space called 
variable dimension space whose dimension changes. This new kind of space 
brings the notion of space and the notion of dimension together. 
The paper is organized as follows. In section [21 a very general adaptive system 
concept called meta- dynamical adaptive system is given by using an extension 
of Kalman systems. Section [3] is dedicated to the Douady's Rabbit fractal 
implemented as a meta-dynamical adaptive system. In section HJ the concept 
of variable dimension space leads to a special model of adaptive system: the 
adaptive differential equation. Finally, we use our concepts of adaptive systems 
to describe a biological example of a defense-attack model in section [51 



2 Meta-dynamical adaptive system 

It is possible to give a general definition of adaptive systems by using the con- 
cept of meta-dynamical adaptive system developed by one of the author during 
his Ph.D. ( [2P] ). But, in order to be able to present our formal approach, we 
will have a look at the main ideas behind the formalization of the concept of 
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dynamical system by Kalman in [12]. The aim of Kalman's approach is to show 
that some very common mathematical structures plus a few axioms can pro- 
vide a very general framework where the notion of dynamical systems (of all 
kinds) is captured. Time is modelled as an ordered subset of the reals (to cover 
both the continuous and discrete paradigms). The important objects are the 
state set (the variables characterizing the system) and the transition function. 
The transition function defines the trajectory in the state set starting from an 
initial state. Only a few axioms are required to characterize these objects and 
allow them to form a "dynamical system". Among the more important ones 
are direction of time, consistency, composition property and causality. Let us 
recall the fundamental definition of Kalman (see [12]). 

A dynamical system respecting Kalman axioms is defined by the following 
axioms: 

(1) There is a given time set T, a state set X, a set of input values U, a set 
of acceptable input functions Q = {u : T ^ U}, a. set of output values 
Y, and a set of output functions F = {7 : T — Y}, 

(2) (Direction of time) T is an ordered subset of the real|l|, 

(3) The input space Q satisfies the following conditions: 

(a) (Nontriviality) Q is nonempty, 

(b) (Concatenation of inputs) An input segment uj^tiM] is 1^ G restricted 
to ]ti, ^2] n T. If tu, cj' G n and ti < ^2 < t^, there is an uj" G Q such 
that a;|;^_,^] = uj]t„t,] and u;|;^,3] = c^j.^^.^], 

(4) There is given a state-transition function 

(p:TxTxXxn^X. 

whose value is the state x (t) = ip {t, r, Xr, u) E X resulting at time t E T 
from the initial state Xr = x (t) G X at initial time t E T under the 
action of the input u E fl. if has the following properties: 

(a) (Direction of time) is defined for alH > r but not necessarily for 
all t < T. 

(b) (Consistency) (p (t, t, x,uj) = x for all t G T, all a; G X and all u E Q. 

(c) (Composition property) For any ti < t2 < t^, we have 

(^3, ti,X,u) = (p {ts, t2, (p {t2, ti, X, U) , Uj) 

for all X G X and all a; G fi. 

(d) (Causality) If cu, cu' G f2 and uj^r,t] = ^]T,t] then 

ip {t, r, x,uj) = ip (t, T, X, oj') . 

(5) There is given a readout map rj : T x X Y which defines the output 
y (t) = r] {t, X (t)). The map ]r, t] ^ Y given by a t-^ r/ (a, (a, r, x, u)), 

^ Such a general definition could include exotic sets such as fractal cantor sets, for 
instance, in practice, the sets used are part of M or N. 
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cr G ]r, t] is an output segment, that is, the restriction 7],- 1] of some 7 G F 
to]r,t]. 

A dynamical system is referred as S = {T, X, U, Q, Y, T, ip, r]}. 

Now, inspired by the definition of dynamical systems of Kalman, we propose 
a formalization of our meta-dynamical adaptive system. 

Definition 1 A meta-dynamical adaptive system M. is a composite mathe- 
matical concept defined by the following axioms: 

(1) The dynamical level: the suitably indexed set 

is a dynamical system respecting Kalman axioms for all {i,j) & I x J 
where are transition functions on state set Xj/ 

ifij -.ToxToX Xixn^ Xi. 

(2) The meta-dynamical level: let X = \J Xj be the set of all the possible 

states of the system and D = {fi.j}{i,j)eixj be the set of all possible 
transition functions, then there exists a meta-dynamical time Ti and a 
meta-dynamical rn/o 

<l>:TixXxL)->Xxi:) 

such that: 

(a) (Temporal hierarchy) Ti C Tq, 

(b) (Coherence states /transitions) // Xi, v^ii jj = (a^2, V^j2,j2) fort G 
Ti, then Xi G Xj^ ^ X2 G Xi^. 

If t E Ti and = {x,ip), $ is said to be mute at (t,x,(p) else 

{t,x,(f) is a commutation point. 

(3) Evolution rule between dynamical and meta-dynamical levels: 
let Xt = ip{t,ti,xi,uj) with uj an input function, there exists a "meta"- 
transition function 

m-.TixTixXxDxVL^XxD 

such that: 

(a) (Dynamical phase) If ^ is mute on (t,Xt,ip) for all t G [ti,t2[nTi, 
then \1/ is defined between ti and t2 and 

^ It is possible to consider input values at the metadynamical level, then $ is defined 
as follows 

<^:TixXxDxni^XxD 
where ^li = {to : Ti ^ Ui} is the set of acceptable "meta" -input functions. 
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(i) if t2 ^ Ti, then 
(a) else t2 G Ti and 

(b) (Concatenation rule) If there exists t2 G]ti,t3[ such that \l/ is defined 
between ti and t2 and between t2 and t^ then \l/ is defined between ti 
and ^3 and 

^(t3,ti,a;i,v?,cj) = '^{t3,t2,'^it2,ti,xi,ip,uj),uj). 

(c) (Stopping rule) \E' is defined between ti and t2 if, respecting previous 
axioms, there is only a finite number of commutation points in [ti, ^2] ■ 

The meta-dynamical rule in point 2) can operate at instants for which the 
system is defined and not necessarily at all of them (see point 2. (a)); actually 
a higher level is usually slower. Moreover, as the meta-dynamical rule can 
change the state and the transition function, we have to consider that both 
match well: the resulting state has to belong to the state set on which the 
new transition function operates (point 2.(6)). In point 3. we describe how 
dynamics and the meta-dynamical rule combine together to make the system 
change with time. 

To take up notions mainly used in social science, our dynamical rule can be 
seen as a heterarchical level and our meta-dynamical level as a hierarchical 
level. Let us recall that a heterarchy is a network of elements which share the 
same "horizontal" position level in a hierarchy. Each level in a hierarchical 
system is composed of a heterarchy which contains its constituent elements 
(see [HI). 

Remark 2 In the second case of (3. a), when t2 G Ti, we deliberately con- 
sider that {t2,^ {xt2,^)) is mute. We do not consider the case where the 
meta-dynamics would "rebound" and have several commutations at the same 
time. If the system has several commutations, it is always possible to con- 
sider this set of commutation as one, with the final state of the last com- 
mutation (provided that, we know from 3.(c) that the number of commuta- 
tions is finite). Axiom 3.(c) is to avoid Zeno-style system with an infinite 
number of commutations in a finite amount of time (e.g. with a quantity 
x{t) = sin (]— ^) with a commutation each time x{t) = on [1 — e, 1 + e]). 

Remark 3 The index set J is associated with the change of dynamics in the 
same state space. We recognize here the framework of hybrid systems (see 
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Remark 4 It is also important to differentiate between the continuous or dis- 
crete dynamics and the continuous or discrete meta- dynamics because the con- 
fusion is easy. The first case is well known. The second one shows the difference 
between a meta- dynamical time Ti which is continuous (see section\^ and a 
meta- dynamical time which is discrete (see section\B^. 



Fig. 1. Illustration of some axioms of a metadynamical system 

Example 5 On Fig. [T], we can see a meta- dynamical system in action. At 
= tc, we pass from a 2— dimensional state set to a 3— dimensional one. 
The dynamics are continuous, so Tq = [to, tend] o-nd the meta- dynamical time 
set Ti is discrete: Ti = {ti,t2,t3, . . . ,ts} . Since the meta-dynamics is mute 
on ti, t2, ts and t^, the evolution from to to t^ is a purely dynamical phase, 
ended by a commutation (Axiom 3. (a) ). Evolution from t^ to tend is also purely 
dynamical (the meta-dynamics is mute on t^, tr andtg. The junction between 
the two of them is made by using the concatenation rule. As there is only one 
commutation, the system is defined. 

Now, it is possible to expand some classical notions of dynamical system, as 
the notion of trajectory. 

Definition 6 Let M. be a meta- dynamical adaptive system with the meta- 
transition function ^ , if for t eTi, 

"^{t, to, xo, ipo, uj) = ipt) 
with Xt G Xt then {Xt,(pt) is a meta-state of Ai in t. 

The sequence {Xt^,ipti)o<i<no<n<+oo of meta-state of M. such that for all t G 
[ti,ti+i[, {Xt,ipt) = {Xt^,ipn) and{Xt.,Lpt^) ^ (Xu+i.^ti+i) is called meta-oihit 
OfM. 

If n < +00, {Xt^,{pt^) is said to be an absorbent meta-state. 

The sequence (^«, -^ti, </'tJo<i<n o<n<+oo satisfying the same conditions is called 
a trajectory of Ai. 

One can extend definition [T] to stochastic systems. 
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Definition 7 Let us consider the dynamical rule of definition^ with the fol- 
lowing meta- dynamical rule 

p$ : Ti X X X D ^ (X X D ^ [0, 1]) 

where is a probability distribution on X xD which represents the probability 

that {x,(f) becomes {xf,!ff) att. The stochastic evolution rule between dynam- 
ical and meta- dynamical levels is given by 

p"^ : Ti X Ti X X X D X n ^ {X X D ^ [0,1]) 

where p'^ is a probability distribution on XxD which represents the probability 

p^ {ti,t2,X,(p,Uj) ■ iXf,(pf) 

that {x,ip) becomes {xf,(pf) between ti and t2 with the input function oj. The 
properties of p^ and p"^ are the same as $ and given m[21 Such a system 
is called a stochastic meta-dynamical adaptive system. 

Tlie interest of definition [7| is its generality. It is for example well adapted to 
model a large number of complex adaptive systems, in particular the biological 
model we develop in section [5l 



3 Algorithm for the Douady's Rabbit fractal 



The goal of this example is to give an algorithm based on definition [Hallowing 
to describe the Douady's Rabbit Fractal. In this example, the set I and thus 
the family {Xjj.gj of definition [T] is not countable. The Riemann Sphere §^ 
is mapped one-to-one onto the extended complex plane Coo = CU {oo} by 
stereographic projection. Let us recall some definitions. Let K = M or C and 
E be the unit ball of K. If / : K ^ IK is a function, x G X is periodic of period 
n G N if /"(a;) = x and for all A; G {1, .., n — 1}, f^{x) ^ x. For x periodic of 
period ra, the cycle 0{x) is the set 



0(x) = {x,/(x),...,r-i(x)} 



and its cardinal is n. Moreover, if / is differentiable at X ^ X IS stable, quasi 



stable or unstable if 



X 



< 1, 



stable, X is attractive if there exists an interval 
for all a; G 



1 or 



X 



> 1. If X is 



V strictly containing x so that 



Fix') 



X. 
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Now, let us recall the definition of the Julia set (see [T]). Let R he a non- 
constant rational function on The Fatou set of R is the maximal open 
subset of on which {-R"}„gj^ is equicontinuous where i?" = R o . . . o R. 
The Julia set Jr of R is the complement of the Fatou set on The filled in 
Julia set of a function R is all the points which are not attracted to the 
super-attracting fixed point at infinity, that is 

Kr = {z eC: {z) 00} . 

This closed set includes the Juha set as its boundary, Jr = OKr. The escape 
set Ir of a function R is all the points that "escape" to infinity, that is 

Ir = {;2 G Coo : i?" (z) ^ 00} . 

If i? is a polynomial of degree 2, Jr is called a quadratic Julia set. The following 
result can be found in [B]: if P is a polynomial of degree d > 2, then Jp is 
closed and dense within itself. 
Here, we are interested in quadratic Julia sets with 

P,iz) =z^ + c (1) 

where c G C. For small values of c, the Julia set is distorted by varying degrees 
from the unit circle, in these cases the Julia set has an infinite length. For 
large values of c, the Julia set becomes an infinite set of totally disconnected 
points, often said to be dust like (in the sense of Cantor). In the quadratic 
case there are only these two possible -connected and disconnected- types. The 
Mandelbrot set is the space containing the value c for which the associated 
Julia set is connected. It is generated by the quadratic sequence 

Zn+l = + C 

(2) 

Zo = C 

\ 

A complex point z = a + ib ^ C will be denoted by (a; 6). The Douady's 
Rabbit fractal is a Mandelbrot set with c = (— | + |(a + 6); ^(a — 6)) where 

a = ^f^^^ and b = so c ~ (-0.12256; 0.744862) (see [7]). As 

the points get closer to the Douady's Rabbit fractal, the speed of convergence 
becomes slower. With the choice of c, {c^ + c) + c = so the origin is an 
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attractive cyclill of period 3 of Pc 



O{0) = {0,c,c^ + c}. 

The boundary points move chaotically. Thus, the idea is to change the points 
by using the speed of convergence as an adaptive value. The principle of our 
algorithm is to have a set of points evolved to the fractal boundary. For this, 
one gives a weight to each point. This weight varies with a dynamics which 
"rewards" the most adapted points (the points which are less attracted by limit 
values) and "penalizes" the least interesting ones. When the most efficient 
points reach a certain threshold, they are allowed to be multiplied in their 
neighborhood. The least efficient points disappear when they reach a minimal 
threshold. This system of thresholds which changes the dynamical structure 
is our meta-dynamics. 

Dynamical level (DL) : with the formalism of definition [H we have Tq = N. 
The index set / is the set of all finite sets of points of C (so i is a set of points 
of C and / is not countable). For all i G I , Xi = is the set of the point 

weights of z. At each given point z E C, we attach a selective value 

/i„ : C ^ M+ 

z ^ rmn{\P^{z)\ , \P^{z) - c\ , \P^{z) - - c\ , ^} 

{0, c, + c} is the attractive cycle of period 3, p^^^-^ is referring to the attrac- 
tion to the infinity seeing that P!^(z) oo is equivalent to -^tttt ~^ 0- 

n— >+oo (^) n^+oo 

One increases the weight of the points which are close to the boundary. In order 
to do this, we organize a "competition" by comparing their mutual slowness 



^ Let K = R or C, / : IK ^ K. be function, x G X is periodic of period n G N if 
/"(x) = X and for all k € {1, ..,n — 1}, f^{x) ^ x. For x periodic of period n, the 
cycle 0{x) is the set 

0{x) = {x,f{x),...,r-\x)} 

and its cardinal is n. Moreover, if / is differentiable at x, x is stable, quasi stable 
or unstable if |(/")'(x)| < 1, |(/")'(x)| = 1 or |(/")'(x)| > 1. If x is stable, x is 
attractive if there exists an interval V strictly containing x so that for all x G F 

r(x') ^ x. 
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of convergence. So, the weights of the points z ^ i are given by 



(t + 1) = {t)+ ^ (-2) with 5 {z, q) 

ooz (0) = 



1 if fit (z) > fit (g) 
if /it (z) = fit (q) 
-1 if yUt (z) < fit (g) 



where t eN. Using the notation of section we have the following transition 
function 

if, {t + 1, t, {UJ, (t)},^,) = {UJ, {t + 

for all t G N. For each state space there is only one associated transition 
function, so the use of the index set J is unnecessary. 

Meta-dynamical level (ML): with the formalism of definition [H we have 
Ti = kN. When the weight of a point reach an upper threshold M > 0, the 
point is allowed to give birth to a new point, randomly in its neighborhood. 
When the weight of a point reaches a lower threshold m < 0, the point is 
removed. This can be modelled by the meta-dynamical rule 

where 

i {n + 1) = {z + ge : g is a. random point in B, z E i (n) and Uz (t) > M} U 
{z E i in) : ojz it) > m} 

such that e > given. The weight of the reproduced point is shared between 
itself Uz (t) and the new neighbouring point Uz+ge (t). The other weights are 
kept equal. Let us sum up the evolution rule: 

• • • ^ "fiin) (t) ^ Viin) + 1) • • • + ^) ^ ^i(n+l) {t + k) ^ . . . 

As the time increases, the selective value becomes more and more accurate, 
i.e. i (n) tends to a set of points belonging to the Julia space Jp^ or to the 
empty set when n +oo. This is an interesting point because the calculations 
are concentrated on the fractal boundary. In a lot of classical algorithms, the 

o 

calculation time is squandered for points in the interior Kp^ (see [IS]). Indeed, 
this meta-dynamical adaptive system produces a cloud of points which gather 
round the Julia space Jp^ called the Douady's Rabbit fractal (see Fig. [2]). 

Though we have developed this algorithm in the setting of Julia sets, the 
same framework can be used to explore many complex frontiers, for example 
other fractal structures where the boundary properties are largely unknown. 
Indeed, this kind of algorithm is really interesting because the calculation is 
concentrated on a specific region. 
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Fig. 2. The Douady's Rabbit fractal 
4 Adaptive differential equations 



First, let us recall the classical notion of Hausdorff fractal dimension (see [17j). 
Let £^ be a subspace of a metric space Ai and p a positive number, one defines 
TZp as the set of all coverings {Bi, pi)^ of S by balls Bi with diameter < pi < p. 
For each positif number a, one denotes: 

H^{S) = znf{Y^^pt:{B,,p,)^en,}. 

H"' (S) = lim if ° (S) is called the a— dimensional Hausdorff measure of £ and 
belongs to [0, +cxd]. Then, let 

dim{£) = m/ {a > : if" {£) = 0} , 

it is the Hausdorff dimension of S. On the one hand, the Hausdorff dimension 
is defined for all metric spaces. On the other hand, in the case of a classical 
space (non fractal), it is identical to its topological dimension (for example the 
Hausdorff dimension of R" is n). In the case of a simple linear fractal, such 
as fractals with internal homothetia obtained by an homothetic iteration with 
constant ratio, the Hausdorff dimension is equal to the homothetic dimension 
dimh (S) given by: 

dimh (S) = = logk{n) 

where n is the number of subsets obtained during the homothetic process of 
reduction with ratio ^ (see [5]). For more information on dimension theory 
the reader may refer to [19] or |10j . 

The variable dimension space is defined as follows: 

Definition 8 Let M. he a metric space and A a parameter space. One defines 
two maps d : A — i> [0, +oo] , A y-^ d{\) and F : A 2^ where F{X) verifies 
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dimF{X) = d{X) and 2-^ is the family of non-empty subsets of A4. F{A) is 
a set of variable dimension spaces. If F{A) is a totally ordered set for the 
inclusion, one calls {F, d) a variable dimension space, d is the dimension 
function and d{F{X)) < d{A4) for all X E A. If A is a topological space then 
F is a set valued function and we may take account of the regularity of d. The 
variable dimension space [F, d) is continuous if d is continuous. Moreover, if 
A is an ordered space then {F,d) is increasing (respectively decreasing^ if d is 
increasing (respectively decreasing). 

To illustrate these definitions, we can consider the following example: 

Example 9 Let | < A < |, one defines four similarities from the family K'^ 
of the compact subset of the square c = [0, 1]^ with value in by 

• Si,A {x) = Xx, 

/l 

V 



S2,X (x) = X 



X 



X + 



A 




S3,A (x) = X 




I 



X + 



and 



Sa\ [x] = Aa; + 



1 - A 




Then, one defines the function Qx '■ K'^ — ^ K'^ , x i— >• si^\ (x)Us2,a (a^)Us3,A (a;)U 
Sax ix). 



Ill, J X-l/4 




0,0 



X,0 



1-X, 



1,0 



Fig. 3. Van Koch snowflake generator 

being contracting, one defines F (A) as the fixed poin{^ of this function 
for the Hausdorff distance. In JJE/, one may find the following result: the 
dimension of F (A) is d (A) = t-tty ■ So, (F, d) is a continuous increasing 



^ To prove this, one uses two classical results: the family of subsets of a Banach space 
endowed with the Hausdorff distance is also a Banach space and if the similarities 
are contracting in K^, then the function J7a is also contracting in K'^ with the 
Hausdorff distance [20l. 
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variable dimension space on A 



1 1 

4' 2 



. A similar construction allows us to 



turn continuously a n— dimension space into a {n + 1) — dimension space. 

An interesting application of the above concept arises when the variable di- 
mension space is the state space of the solutions of a differential equation 
over a period of time. The union of the solutions then "moves" on a variable 
dimension space. We talk about an adaptive differential equation. 

Definition 10 Let X be a Banach space of finite dimension called the possible 
state space. Consider 

(1) a subdivision {tjjjgp^ o/]R_|_, 

(2) an application d : M+ x X ^ N, (t, y) ^— (i (t, y) such that d (0, y) is given, 

(3) an application g:R+x M'^*'^(^) W^^^i^) ^ 

such that one has the system 

x^it) = fi{t,Xi{t)),t G [ti,ti+i[,x,{t) G R''^'-y^\t e N 
Xi (tr) = g {ti, yi) 

where yo = Xq (0), y^ = lim (t) for i G N*, Xi (.) is the right derivative of 
Xi (.) and 

is a family of applications. Such a system is called an adaptive differential 
equation. A trajectory of the system ^ is a family 

x{t) = {xi (t) : t G [ti,ti+i[}.^^. 



One may notice that dim{g (ti,yi)) = d(ti,yi). This concept of adaptive differ- 
ential equation is not a succession of ordinary differential equations because 
the initial condition of each system i depends on the system i — 1. The ordi- 
nary differential equations represent where d is constant. 
The adaptive differential equation is a reductive approach of the meta-dynamical 
adaptive system to a system described by a differential equation whose meta- 
dynamics is discrete and where dimX < -|-cxo. 



5 Application to a biological model 



Here, we want to model the influence of the bacterium-phage interaction on 
the co-evolution of the populations of bacteria and phages. In the following we 
will need the definition of the Hamming distance. Let n G N, the Hamming 
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distance is the function dn '■ {0, 1} 



n 



N defined by 



dn {si,S2) = 




wfiere si, S2 € {0, 1}". dn represents the number of differing bits between the 
two binary strings Si = s\sl . . . and S2 = sls2 . . . 

The attack of a bacteria population by phages is assumed to be done by the 
lysis process: a phage hangs on the surface of a bacterium cell, injects its 
DNA in it and then forces the bacterium to yield its own replicas inside the 
cell. When the cell is full, it bursts, releasing a huge quantity of copies of the 
infecting phage. The efficiency of the attack (i.e. the probability of success 
of the infection), depends on the couple bacterium-phage. One of the other 
characteristics of bacterial and phages populations are their high variabilities. 
They frequently mutate, creating new populations with new properties. Such 
a system has two dynamics to be taken into account: the dynamics of the 
populations of bacteria and phages and the meta-dynamics of evolution geared 
by mutations and extinctions. It is without any doubt a model which is not 
in the scope of the classical theory of dynamical systems, but our formalism 
applies well to it. 

Dynamical level: it is made up of a set of ordinary differential equations. 
Let us consider the following system 



with S the concentration of nutrient, Bi the different bacteria strains with 
< i < Ub and Pj and Ikj the different phages strains with < j < np. 
To each population Bi and Pj is associated a binary string G {0, 1}" and 
G {0, 1}". These binary strings code for the properties of attack (for the 
phages) or defence (for the bacteria) facing the infection (see [2] for the bi- 
ological discussion). The populations are then characterized by the different 
concentrations and the two lists of bit strings. The set L of all the pairs of 
lists is taken as our second index (the one called J in definition [1]). 

The model is described by a modified version of Mosekilde equations (see 
jl5]). This set of ordinary differential equations describes the interactions of 
bacterial populations Bi and phage populations Pj in a chemostat. Bi and Pj 
also symbolize the concentration of these populations and S the concentration 
of the nutrient. The process of infection of bacteria by phages is modelled by 
three infection stages Iij, I2J and I^j. One associates with the ecosystem 
dl]) the state space Xi^nt+iup = IR^"*""''^^"'''. Over a period of time without 
appearance or disappearance of any strain of bacteria or phages, the dynamical 
evolution of the system (jl]) is modelled by the evolution of the concentrations 
of S, the different bacteria Bi and phages Pj, and three infection stages Ikj- 



{S, {-Bi}o< i<nh 1 



{Pj, hj, hj, h,j}o<j<nJ 



(4) 
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With the formahsm of definition [H that means that Tq = M and for a fixed 
/ G L the transition function ipi^nb+inp,i is the integration of the following set 
of differential equations 



= ^ - i?. E au,,P, - pB, 



dl_ 

dt 
dt 
dt 



dPj ^ oPh 
dt T 



i=l 

l{h,3 - hi) - Phi 

' nt rip 3 



dS 
dt 



P, E + E E 

\i=l i=i fc=i ^ 

"6 



wySBj 

K+S 



(5) 



with < z < -n,;,, < J < rip, p the rate of dilution (p = 0.0045 min~^), 
Kand z/ respectively the saturation term and the growth from the Monod 
equation (k = 10pg.ml~^, u = 0.024 min~^), a the theoretical adsorption 
constant depending on phage and bacterium size {a = 10~^ml. min~^), t 
a time constant (r = 30 min), (3 the number of copies of phage j released 
during the burst of the infected bacterial cell (/? = 100), a the continuous 
supply of substrate (a = 10pg.ml~^), 7 the amount of nutrient consumed in 
each cellular division (7 = O.Oln^f) and finally Uij the probability of infection 
of Bi by Pj which depends on the similarity between bit string (attached to 
bacterium population Bi) and (attached to phage population Pj) as follows 

= V — ~) 

with He the size of the binary string. 

Here, we are not interested in the identification of the biological dynamical 
level which can be found in [2]. With the given size of a binary string Uc, 
there exists a finite number of possible dynamical systems, here differential 
equations. Indeed, the number of possible different populations of bacteria is 
equal to the number of parts of the set of binary strings of size Uc- So, there are 
2"^""" possible populations of bacteria. For the same reason, one deduces that 
the possible number of populations of phages is the same and thus the total 
number of possible state spaces is 2"^"" x 2^"'' = 2'^""'^^. So, with the notation of 
definition [H it means that card (J) = 2'^""'^^. If we take for example Uc = 10, 
one has 2^°^^ ~ 3.6 x 10^°^ possible state spaces. Theoretically, one can consider 
a 3.6 X 10'^°^ dimensional space to embed the system ([5]). Nevertheless, when it 
comes to numerical simulations, such a big system is impossible to deal with. 
One thus understands the need of an adaptive system to describe the system 
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Meta-dynamical level (meta-dynamical adaptive system point of 
view): this is the main difference with the model of Mosekilde given in [151 
Chapter A] which is not evolutive. 

Proposition 11 Consider a small interval of time At, the adaptive changes 
of the system ^ are given by the following mechanisms 

A (t)^ 

p^{t + At, t, 6, (p2+n,+4np,h) " (^1' ■■■^^k, (pi+ks+n,+4np,h) = ^'^^^ 

and 

'^{t + At, t, (5, (p2+nt,+4np,h , ^s) = (/?2-E4+nt,+4np,«2 (7) 

where 



(1) formula ^ is the probability that the species S (a bacterium B or a 
phage P) gives birth to k > 1 mutante strains 6i,. . . ,6k on [t,t + At[ 
with A (t) = ^Pe, Pe the probability that a small group of mutant species 
of size 6e gradually replacing the species of the parent population 6. (pij 
is the function defined by equations ^ for j = li,l2 and ks = k if 6 is a 
bacterium and kg = Ak if 6 is a phage, 

(2) formula ^ gives the determinist rule of the extinction of the species S 
which depends on a given threshold rris, es = 1 if 5 is a bacterium and 
es = 4: if 6 is a phage. 

p'^ is the stochastic transition function of definition [7] which governs the meta- 
dynamical rule of the appearance. is the transition function of definition [T] 
which governs the meta-dynamical rule of the extinction. 

Proof. Consider the probability PBiitjk) that the population Bi gives birth 
to k mutante strain on [t, t + At[ and the probability pe that a small group of 
mutante bacteria of size Be is gradually replacing the bacteria of the parent 
population Bi. Such a reasoning gives a binomial probability for 

with n(t) = When n(t) is large, one may approximate the binomial 

probability by the Poisson probability 

p^,(t,A:)^e-^wA^ (8) 

with A (t) = ^^Pe- Suppose that the birth of all the populations to a Ham- 
ming distance of one (only one bit is different) is equiprobable. The birth of a 
mutant strain results in the change of a group of bacteria (which is a part of a 
parent population) of size -Bg, in a population group with new characteristics. 
If this population group already exists, the mutant population is added to 
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it. The same mechanism governs the modelhng of the phage mutation. For 
the extension, there exists a threshold below which the extinction of the 
population is certain. Every population under a given threshold (different for 
bacteria and phages mp) is removed from the system. Thus, the extinction 
is determinist. ■ 

We have defined the macroscopic birth of a mutant population as an event 
occurring on an interval [t,t + At[. There exists a time set Ti where the 
system commutes. This commutation depends on the state of the system 
and on the coefficient Pb pp respectively coefficient of the Poisson law 
of the bacterium and the phage. In order to simplify the model and to make 
it computable, one may suppose that 

Ti = to + i^t 

even if in practice, the mutations have no reasons to be periodically defined. 
This allows to define, with the previous notations, the meta-dynamical rule 

p^{to + iAt, 6, V?2+n6+4np,h) " {^U ■■■,6k, ipi+ks+nt+4np,l2) = ^'^^^^ 

Then, at each step At there are four possible commutations: 

• birth of a new bacterial strain: a variable -Bnj,+i is added, the dimension of 
the system ([5]) increases by one, 

• birth of a new phagical strain: four variables are added: Pnp+i, h,np+i, h,np+i 
and h^np+i, the dimension of the system ([5]) increases by four, 

• extinction of a bacterial strain: the concerned variable is removed, the di- 
mension of the system ([5]) decreases by one, 

• extinction of phagical strain: variables of the concerned phage are removed, 
the dimension of the system ([5]) decreases by four. 

The different possible state spaces resulting from a commutation are given on 
figure m 

If more than one event occurs at each moment t E Ti, one composes the pos- 
sible change of state spaces (for example, the extinction of a phage combined 
with the birth of a bacterial strain decreases the dimension of the system ([5]) 
by 4- 1 = 3). 

Meta-dynamical level (adaptive differential equation point of view) : 

we have described the biological model as a stochastic meta-dynamical adap- 
tive system by using definition [71 This modelling corresponds to a stochastic 
view of the system (j4]) where the space of all stochastic realizations is infi- 
nite. There exists a family of points {tjjjgp^ where the system ([5]) commutes. 
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Appearance of a new 
strain of bacteria 




Fig. 4. Evolution of the state space with Ub and Up respectively the number of 
bacterial and phage strains before the transition. 

We have supposed that this commutation is given by ti = to + i^t. Then, 
at each step At the system may commute. To see a possible evolution of the 
system we choose a stochastic realization at each point t, (only for the ap- 
pearance of a strain because the extinction of a strain is determinist) g{ti,yi) 
where = lim Xi-i (t) with Xi-i(t) the solution of the system ([5]) on [tj_i,ti[. 



It describes the determinist evolution of the system 
a stochastic realization is a function 



By extension, for us 
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X 



ndim{X) 



ndim(X) 



defined at least on {iti,yi)}-^^. For a given realization g, the system (jlj) may 
be modelled by an adaptive differential equation whose equations are given 
by the system ([5]). This modelling belongs to the variable dimensional space 
j^ni,+4np+i ^Yiere + Aup + 1 depends on g and follows the rule given on figure 
m A detailed study of this biological model with implementation can be found 
in [3]. 

In this example, one sees that the framework of the adaptive differential equa- 
tions corresponds to the case of a transition function defined by a set of dif- 
ferential equations and a meta-dynamical time reduced to a discrete subset of 



6 Conclusion 



We started our research with a biological system whose dynamics changes 
in different dimensions. Seeing that there was no mathematical framework 
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to describe such a system, we have developed a mathematical tool following 
Kalman's dynamical system called meta-dynamical adaptive system which 
was appropriate to give a constructive algorithm for some fractals. It was 
also adapted to describe and analyze our biological system. However, we have 
found our tool too general and we decide to develop a special tool for dif- 
ferential equations. The new system called "adaptive differential equation" 
is not a succession of differential equations because the initial condition of 
each system depends on the previous system and gives the new dimension 
of the following system. This model allows to describe a system of changing 
differential equations, in particular the stochastic realization of a stochastic 
meta-dynamical adaptive system. The last tool we use is the "variable dimen- 
sion space" . This new kind of space links the notion of space and dimension in 
a changing dynamics. We think that our work will contribute to understand 
the huge number of complex systems where the espace of exploration is too 
big to be investigate with classical means. 
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1 Introduction 



In the two last decades, there has been much interest in the study and for- 
malization of complex adaptive systems (see [?] and [?]). Many different ap- 
proaches have been proposed: artificial chemistries, evolutionary formalism, 
cellular automata. However, from a more theoretical viewpoint, few mathe- 
matical formalisms exist for adaptive systems. Though we may cite the chap- 
ter "Categorical System Theory" proposed by A.H. Louie in [?] who discuss 
the relationship between natural and formal systems, most of the studies are 
mainly simulation-based (see [?] and [?]). With the recent evolution of physics 
and biology, a general mathematical formalism for adaptive systems would be 
very useful. Actually, behind their apparent heterogeneity, most of the adap- 
tive systems share one important feature: they are dynamical objects whose 
structures are sometimes modified by a top level automation-like rule. On the 
basis of this observation, we have built a two-level formalism that helps us to 
design an algorithm for fractals and a biological model of co-evolution in a 
bacterium-phage system. 

The main goal of this article is to present a concept of adaptive systems, gen- 
eral enough to be used in different fields (mathematics, physics or biology) and 
to apply this framework to a fractal implementation and a biological system. 
Besides, a new concept of spaces is developed. Since the second part of the 
19*^ century, various definitions of dimensions, in particular those developed 
by Cantor and Peano, have appeared (see [?]). These new definitions lead to 
the concept of a fractal model and with it, a global view of spaces (see [?]). In 
order to study adaptive systems, we will use a sort of "adaptive" space called 
variable dimension space whose dimension changes. This new kind of space 
brings the notion of space and the notion of dimension together. 
The paper is organized as follows. In section [21 a very general adaptive system 
concept called meta- dynamical adaptive system is given by using an extension 
of Kalman systems. Section [3] is dedicated to the Douady's Rabbit fractal 
implemented as a meta-dynamical adaptive system. In section HJ the concept 
of variable dimension space leads to a special model of adaptive system: the 
adaptive differential equation. Finally, we use our concepts of adaptive systems 
to describe a biological example of a defense-attack model in section [51 



2 Meta-dynamical adaptive system 

It is possible to give a general definition of adaptive systems by using the con- 
cept of meta-dynamical adaptive system developed by one of the author during 
his Ph.D. ([?,?]). But, in order to be able to present our formal approach, we 
will have a look at the main ideas behind the formalization of the concept of 
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dynamical system by Kalman in [?] . The aim of Kalman's approach is to show 
that some very common mathematical structures plus a few axioms can pro- 
vide a very general framework where the notion of dynamical systems (of all 
kinds) is captured. Time is modelled as an ordered subset of the reals (to cover 
both the continuous and discrete paradigms). The important objects are the 
state set (the variables characterizing the system) and the transition function. 
The transition function defines the trajectory in the state set starting from an 
initial state. Only a few axioms are required to characterize these objects and 
allow them to form a "dynamical system". Among the more important ones 
are direction of time, consistency, composition property and causality. Let us 
recall the fundamental definition of Kalman (see [?]). 

A dynamical system respecting Kalman axioms is defined by the following 
axioms: 

(1) There is a given time set T, a state set X, a set of input values U, a set 
of acceptable input functions Q = {u : T ^ U}, a. set of output values 
Y, and a set of output functions F = {7 : T — Y}, 

(2) (Direction of time) T is an ordered subset of the real|l|, 

(3) The input space Q satisfies the following conditions: 

(a) (Nontriviality) Q is nonempty, 

(b) (Concatenation of inputs) An input segment uj^tiM] is 1^ G restricted 
to ]ti, ^2] n T. If uj,Lj' E Q and ti < ^2 < ^3, there is an u" G Q such 
that a;|;^_,^] = uj]t„t,] and u;|;^,3] = c^j.^^.^], 

(4) There is given a state-transition function 

(f-.TxTxXxn^X. 

whose value is the state x (t) = ip {t, r, Xr, u) E X resulting at time t E T 
from the initial state Xr = x (t) E X at initial time t E T under the 
action of the input u E fl. if has the following properties: 

(a) (Direction of time) is defined for alH > r but not necessarily for 
all t < T. 

(b) (Consistency) (p (t, t, x,uj) = x for all t E T, all x E X and all u E Q. 

(c) (Composition property) For any ti < t2 < t^, we have 

(^3, ti,X,u) = (p {ts, t2, (p {t2, ti, X, U) , Uj) 

for all X G X and all u E Q. 

(d) (Causality) If u!,u!' E fl and uj^r,t] = ^]T,t] then 

ip {t, r, x,uj) = ip (t, T, X, oj') . 

(5) There is given a readout map rj : T x X Y which defines the output 
y (t) = r] {t, X (t)). The map ]r, t] ^ Y given by a t-^ r/ (a, (a, r, x, u)), 

^ Such a general definition could include exotic sets such as fractal cantor sets, for 
instance, in practice, the sets used are part of M or N. 
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cr G ]r, t] is an output segment, that is, the restriction 7],- 1] of some 7 G F 
to]r,t]. 

A dynamical system is referred as S = {T, X, U, Q, Y, T, ip, r]}. 

Now, inspired by the definition of dynamical systems of Kalman, we propose 
a formalization of our meta-dynamical adaptive system. 

Definition 1 A meta-dynamical adaptive system M. is a composite mathe- 
matical concept defined by the following axioms: 

(1) The dynamical level: the suitably indexed set 

is a dynamical system respecting Kalman axioms for all {i,j) & I x J 
where are transition functions on state set Xj/ 

ifij -.ToxToX Xixn^ Xi. 

(2) The meta-dynamical level: let X = \J Xj be the set of all the possible 

states of the system and D = {fi.j}{i,j)eixj be the set of all possible 
transition functions, then there exists a meta-dynamical time Ti and a 
meta-dynamical rn/o 

<l>:TixXxL)->Xxi:) 

such that: 

(a) (Temporal hierarchy) Ti C Tq, 

(b) (Coherence states /transitions) // Xi, v^ii jj = (a^2, V^j2,j2) fort G 
Ti, then Xi G Xj^ ^ X2 G Xi^. 

If t E Ti and = {x,ip), $ is said to be mute at (t,x,(p) else 

{t,x,(f) is a commutation point. 

(3) Evolution rule between dynamical and meta-dynamical levels: 
let Xt = ip{t,ti,xi,uj) with uj an input function, there exists a "meta"- 
transition function 

m-.TixTixXxDxVL^XxD 

such that: 

(a) (Dynamical phase) If ^ is mute on (t,Xt,ip) for all t G [ti,t2[nTi, 
then \1/ is defined between ti and t2 and 

^ It is possible to consider input values at the metadynamical level, then $ is defined 
as follows 

<^:TixXxDxni^XxD 
where ^li = {to : Ti ^ Ui} is the set of acceptable "meta" -input functions. 
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(i) ift2 ^ Ti, then 

^'(^2,^1, 2:1,93,0;) = (a;t2,<^), 

(ii) else t2 G Ti and 

^'(^2,^1, 2:1,93, a;) = ^{xt^,(p). 

(b) (Concatenation rule) If there exists t2 ^]ti,t-i[ such that ^ is defined 
between ti and t^ and between t2 and t^ then ^ is defined between ti 
and ts and 

'i>{t3,ti,Xi,(fi,Uj) = '^{t3,t2,'^{t2,ti,Xi,(fi,Uj),Uj). 

(c) (Stopping rule) \t' is defined between ti and t2 if, respecting previous 
axioms, there is only a finite number of commutation points in [^1,^2]- 

The meta-dynamical rule in point 2) can operate at instants for which the 
system is defined and not necessarily at all of them (see point 2. (a)); actually 
a higher level is usually slower. Moreover, as the meta-dynamical rule can 
change the state and the transition function, we have to consider that both 
match well: the resulting state has to belong to the state set on which the 
new transition function operates (point 2.(6)). In point 3. we describe how 
dynamics and the meta-dynamical rule combine together to make the system 
change with time. 

To take up notions mainly used in social science, our dynamical rule can be 
seen as a heterarchical level and our meta-dynamical level as a hierarchical 
level. Let us recall that a heterarchy is a network of elements which share the 
same "horizontal" position level in a hierarchy. Each level in a hierarchical 
system is composed of a heterarchy which contains its constituent elements 
(see [?]). 

Remark 2 In the second case of (3. a), when t2 G Ti, we deliberately con- 
sider that (t2,^ {xt2,^)) is mute. We do not consider the case where the 
meta-dynamics would "rebound" and have several commutations at the same 
time. If the system has several commutations, it is always possible to con- 
sider this set of commutation as one, with the final state of the last com- 
mutation (provided that, we know from 3.(c) that the number of commuta- 
tions is finite). Axiom 3.(c) is to avoid Zeno-style system with an infinite 
number of commutations in a finite amount of time (e.g. with a quantity 
x{t) — sin (yjzi) '^'^th a commutation each time x(t) — on [1 — e, 1 -|- e\). 

Remcirk 3 The index set J is associated with the change of dynamics in the 
same state space. We recognize here the framework of hybrid systems (see [?]). 
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Remark 4 It is also important to differentiate between the continuous or dis- 
crete dynamics and the continuous or discrete meta- dynamics because the con- 
fusion is easy. The first case is well known. The second one shows the difference 
between a meta- dynamical time Ti which is continuous (see section\^ and a 
meta- dynamical time which is discrete (see section\B^. 



Fig. 1. Illustration of some axioms of a metadynamical system 

Example 5 On Fig. [T], we can see a meta- dynamical system in action. At 
= tc, we pass from a 2— dimensional state set to a 3— dimensional one. 
The dynamics are continuous, so Tq = [to, tend] o-nd the meta- dynamical time 
set Ti is discrete: Ti = {ti,t2,t3, . . . ,ts} . Since the meta-dynamics is mute 
on ti, t2, ts and t^, the evolution from to to t^ is a purely dynamical phase, 
ended by a commutation (Axiom 3. (a) ). Evolution from t^ to tend is also purely 
dynamical (the meta-dynamics is mute on t^, tr andtg. The junction between 
the two of them is made by using the concatenation rule. As there is only one 
commutation, the system is defined. 

Now, it is possible to expand some classical notions of dynamical system, as 
the notion of trajectory. 

Definition 6 Let M. be a meta- dynamical adaptive system with the meta- 
transition function ^ , if for t eTi, 

"^{t, to, xo, ipo, uj) = ipt) 
with Xt G Xt then {Xt,(pt) is a meta-state of Ai in t. 

The sequence {Xt^,ipti)o<i<no<n<+oo of meta-state of M. such that for all t G 
[ti,ti+i[, {Xt,ipt) = {Xt^,ipn) and{Xt.,Lpt^) ^ (Xu+i.^ti+i) is called meta-oihit 
OfM. 

If n < +00, {Xt^,{pt^) is said to be an absorbent meta-state. 

The sequence (^«, -^ti, </'tJo<i<n o<n<+oo satisfying the same conditions is called 
a trajectory of Ai. 

One can extend definition [T] to stochastic systems. 
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Definition 7 Let us consider the dynamical rule of definition^ with the fol- 
lowing meta- dynamical rule 

p$ : Ti X X X D ^ (X X D ^ [0, 1]) 

where is a probability distribution on X xD which represents the probability 

that {x,(f) becomes {xf,!ff) att. The stochastic evolution rule between dynam- 
ical and meta- dynamical levels is given by 

p"^ : Ti X Ti X X X D X n ^ {X X D ^ [0,1]) 

where p'^ is a probability distribution on XxD which represents the probability 

p^ {ti,t2,X,(p,Uj) ■ iXf,(pf) 

that {x,ip) becomes {xf,(pf) between ti and t2 with the input function oj. The 
properties of p^ and p"^ are the same as $ and given m[21 Such a system 
is called a stochastic meta-dynamical adaptive system. 

Tlie interest of definition [7| is its generality. It is for example well adapted to 
model a large number of complex adaptive systems, in particular the biological 
model we develop in section [5l 



3 Algorithm for the Douady's Rabbit fractal 



The goal of this example is to give an algorithm based on definition [Hallowing 
to describe the Douady's Rabbit Fractal. In this example, the set I and thus 
the family {Xjj.gj of definition [T] is not countable. The Riemann Sphere §^ 
is mapped one-to-one onto the extended complex plane Coo = CU {oo} by 
stereographic projection. Let us recall some definitions. Let K = M or C and 
E be the unit ball of K. If / : K ^ IK is a function, x G X is periodic of period 
n G N if /"(a;) = x and for all A; G {1, .., n — 1}, f^{x) ^ x. For x periodic of 
period ra, the cycle 0{x) is the set 



0(x) = {x,/(x),...,r-i(x)} 



and its cardinal is n. Moreover, if / is differentiable at X ^ X IS stable, quasi 



stable or unstable if 



X 



< 1, 



stable, X is attractive if there exists an interval 
for all a; G 



1 or 



X 



> 1. If X is 



V strictly containing x so that 



Fix') 



X. 
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Now, let us recall the definition of the Julia set (see [?]). Let R he a, non- 
constant rational function on S^. The Fatou set of R is the maximal open 
subset of on which is equicontinuous where R^ = R o . . . o R. 

The Julia set Jr of R is the complement of the Fatou set on S^. The filled in 
Julia set Kr of a function R is all the points which are not attracted to the 
super-attracting fixed point at infinity, that is 

KR^izeC-.R"" (z) 00} . 

This closed set includes the Julia set as its boundary, Jr = OKr. The escape 
set Ir of a function R is all the points that "escape" to infinity, that is 

lR = {zeC^: R"" {z) ^ 00} . 

If i? is a polynomial of degree 2, Jr is called a quadratic Julia set. The following 
result can be found in [?]: if P is a polynomial of degree d > 2, then Jp is 
closed and dense within itself. 
Here, we are interested in quadratic Julia sets with 

P,{z) = z^ + c (1) 

where c G C. For small values of c, the Julia set is distorted by varying degrees 
from the unit circle, in these cases the Julia set has an infinite length. For 
large values of c, the Julia set becomes an infinite set of totally disconnected 
points, often said to be dust like (in the sense of Cantor). In the quadratic 
case there are only these two possible -connected and disconnected- types. The 
Mandelbrot set is the space containing the value c for which the associated 
Julia set is connected. It is generated by the quadratic sequence 



Zn+l — Z^ + C 

< 

Zo^C 

A complex point z = a + ib & C will be denoted by (a; 6). The Douady's 
Rabbit fractal is a Mandelbrot set with c = (— | + |(a + 6); ^(a — 6)) where 

a = ^ ^^+^ and b = so c ~ (-0.12256; 0.744862) (see [?]). As 

the points get closer to the Douady's Rabbit fractal, the speed of convergence 
becomes slower. With the choice of c, (c^ -|- c)^ -|- c = so the origin is an 
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(2) 



attractive cyclill of period 3 of Pc 



O{0) = {0,c,c^ + c}. 

The boundary points move chaotically. Thus, the idea is to change the points 
by using the speed of convergence as an adaptive value. The principle of our 
algorithm is to have a set of points evolved to the fractal boundary. For this, 
one gives a weight to each point. This weight varies with a dynamics which 
"rewards" the most adapted points (the points which are less attracted by limit 
values) and "penalizes" the least interesting ones. When the most efficient 
points reach a certain threshold, they are allowed to be multiplied in their 
neighborhood. The least efficient points disappear when they reach a minimal 
threshold. This system of thresholds which changes the dynamical structure 
is our meta-dynamics. 

Dynamical level (DL) : with the formalism of definition [H we have Tq = N. 
The index set / is the set of all finite sets of points of C (so i is a set of points 
of C and / is not countable). For all i G I , Xi = is the set of the point 

weights of z. At each given point z E C, we attach a selective value 

/i„ : C ^ M+ 

z ^ rmn{\P^{z)\ , \P^{z) - c\ , \P^{z) - - c\ , ^} 

{0, c, + c} is the attractive cycle of period 3, p^^^-^ is referring to the attrac- 
tion to the infinity seeing that P!^(z) oo is equivalent to -^tttt ~^ 0- 

n— >+oo (^) n^+oo 

One increases the weight of the points which are close to the boundary. In order 
to do this, we organize a "competition" by comparing their mutual slowness 



^ Let K = R or C, / : IK ^ K. be function, x G X is periodic of period n G N if 
/"(x) = X and for all k € {1, ..,n — 1}, f^{x) ^ x. For x periodic of period n, the 
cycle 0{x) is the set 

0{x) = {x,f{x),...,r-\x)} 

and its cardinal is n. Moreover, if / is differentiable at x, x is stable, quasi stable 
or unstable if |(/")'(x)| < 1, |(/")'(x)| = 1 or |(/")'(x)| > 1. If x is stable, x is 
attractive if there exists an interval V strictly containing x so that for all x G F 

r(x') ^ x. 
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of convergence. So, the weights of the points z ^ i are given by 



(t + 1) = {t)+ ^ (-2) with 5 {z, q) 

ooz (0) = 



1 if fit (z) > fit (g) 
if /it (z) = fit (q) 
-1 if yUt (z) < fit (g) 



where t eN. Using the notation of section we have the following transition 
function 

if, {t + 1, t, {UJ, (t)},^,) = {UJ, {t + 

for all t G N. For each state space there is only one associated transition 
function, so the use of the index set J is unnecessary. 

Meta-dynamical level (ML): with the formalism of definition [H we have 
Ti = kN. When the weight of a point reach an upper threshold M > 0, the 
point is allowed to give birth to a new point, randomly in its neighborhood. 
When the weight of a point reaches a lower threshold m < 0, the point is 
removed. This can be modelled by the meta-dynamical rule 

where 

i {n + 1) = {z + ge : g is a. random point in B, z E i (n) and Uz (t) > M} U 
{z E i in) : ojz it) > m} 

such that e > given. The weight of the reproduced point is shared between 
itself Uz (t) and the new neighbouring point Uz+ge (t). The other weights are 
kept equal. Let us sum up the evolution rule: 

• • • ^ "fiin) (t) ^ Viin) + 1) • • • + ^) ^ ^i(n+l) {t + k) ^ . . . 

As the time increases, the selective value becomes more and more accurate, 
i.e. i (n) tends to a set of points belonging to the Julia space Jp^ or to the 
empty set when n +oo. This is an interesting point because the calculations 
are concentrated on the fractal boundary. In a lot of classical algorithms, the 

o 

calculation time is squandered for points in the interior Kp^ (see [?]). Indeed, 
this meta-dynamical adaptive system produces a cloud of points which gather 
round the Julia space Jp^ called the Douady's Rabbit fractal (see Fig. [2]). 

Though we have developed this algorithm in the setting of Julia sets, the 
same framework can be used to explore many complex frontiers, for example 
other fractal structures where the boundary properties are largely unknown. 
Indeed, this kind of algorithm is really interesting because the calculation is 
concentrated on a specific region. 
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Fig. 2. The Douady's Rabbit fractal 
4 Adaptive differential equations 

First, let us recall the classical notion of Hausdorff fractal dimension (see [?]). 
Let £ be a subspace of a metric space A4 and p a positive number, one defines 
Tip as the set of all coverings {Bi, pi),- of S by balls Bi with diameter < pi < p. 
For each positif number a, one denotes: 

H^{£)^inf{Y^^pt:iB„p,),enp}. 

i£) = hm if" (£) is called the a— dimensional Hausdorff measure of S and 
belongs to [0, +oo]. Then, let 

dim(£) = in/ {a > : (S) = 0} , 

it is the Hausdorff' dimension of £. On the one hand, the Hausdorff dimension 
is defined for all metric spaces. On the other hand, in the case of a classical 
space (non fractal), it is identical to its topological dimension (for example the 
Hausdorff dimension of M"^ is n). In the case of a simple linear fractal, such 
as fractals with internal homothetia obtained by an homothetic iteration with 
constant ratio, the Hausdorff dimension is equal to the homothetic dimension 
dinih (S) given by: 

diruh (S) = 1^1^ = logk{n) 

where n is the number of subsets obtained during the homothetic process of 
reduction with ratio ^ (see [?]). For more information on dimension theory 
the reader may refer to [?] or [?]. 

The variable dimension space is defined as follows: 

Definition 8 Let M. he a metric space and A a parameter space. One defines 
two maps d : [0, +oo] , A ^ d{X) and F : A ^ 2^ where F{X) verifies 
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dimF{X) = d{X) and 2-^ is the family of non-empty subsets of A4. F{A) is 
a set of variable dimension spaces. If F{A) is a totally ordered set for the 
inclusion, one calls {F, d) a variable dimension space, d is the dimension 
function and d{F{X)) < d{A4) for all X E A. If A is a topological space then 
F is a set valued function and we may take account of the regularity of d. The 
variable dimension space [F, d) is continuous if d is continuous. Moreover, if 
A is an ordered space then {F,d) is increasing (respectively decreasing^ if d is 
increasing (respectively decreasing). 

To illustrate these definitions, we can consider the following example: 

Example 9 Let | < A < |, one defines four similarities from the family K'^ 
of the compact subset of the square c = [0, 1]^ with value in by 

• Si,A {x) = Xx, 

/l 

V 



S2,X (x) = X 



X 



X + 



A 




S3,A (x) = X 




I 



X + 



and 



Sa\ [x] = Aa; + 



1 - A 




Then, one defines the function Qx '■ K'^ — ^ K'^ , x i— >• si^\ (x)Us2,a (a^)Us3,A (a;)U 
Sax ix). 



Ill, J X-l/4 




0,0 



X,0 



1-X, 



1,0 



Fig. 3. Van Koch snowflake generator 

ilx being contracting, one defines F (A) as the fixed poin{^ of this function for 
the Hausdorff distance. In [ ?], one may find the following result: the dimen- 
sion of F (A) is d (A) = T-TTY ■ So, (F, d) is a continuous increasing variable 



^ To prove this, one uses two classical results: the family of subsets of a Banach space 
endowed with the Hausdorff distance is also a Banach space and if the similarities 
are contracting in K^, then the function J7a is also contracting in K'^ with the 
Hausdorff distance [?1. 
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dimension space on A = . A similar construction allows us to turn con- 
tinuously a n— dimension space into a {n + 1)~ dimension space. 

An interesting application of the above concept arises when the variable di- 
mension space is the state space of the solutions of a differential equation 
over a period of time. The union of the solutions then "moves" on a variable 
dimension space. We talk about an adaptive differential equation. 

Definition 10 Let X be a Banach space of finite dimension called the possible 
state space. Consider 

(1) a subdivision {tjj.gj^ o/]R+, 

(2) an application d : M+ x X ^ fi, {t,y) ^ d (t, y) such that d (0, y) is given, 

(3) an application g:R+x M*™(^) W^^^i^) ^ 

such that one has the system 

Xiif) = fi{t,Xi{t)),te [ti,ti+i[,xi{t) e R''^'-y^\t e N 
Xi (tr) = g {ti, yi) 

where yo = Xq (0), yi = lim (t) for i G N*, Xi (.) is the right derivative of 
Xi (.) and 

is a family of applications. Such a system is called an adaptive differential 
equation. A trajectory of the system ^ is a family 

x{t) = {xi (t) -.te 



One may notice that dim{g (ti,yi)) = d{ti,yi). This concept of adaptive differ- 
ential equation is not a succession of ordinary differential equations because 
the initial condition of each system i depends on the system i — 1. The ordi- 
nary differential equations represent where d is constant. 
The adaptive differential equation is a reductive approach of the meta-dynamical 
adaptive system to a system described by a differential equation whose meta- 
dynamics is discrete and where dimX < -|-cxo. 



5 Application to a biological model 



Here, we want to model the influence of the bacterium-phage interaction on 
the co-evolution of the populations of bacteria and phages. In the following we 
will need the definition of the Hamming distance. Let n G N, the Hamming 
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distance is the function dn '■ {0, 1} 



n 



N defined by 



dn {si,S2) = 




wfiere si, S2 € {0, 1}". dn represents the number of differing bits between the 
two binary strings Si = s\sl . . . and S2 = sls2 . . . 

The attack of a bacteria population by phages is assumed to be done by the 
lysis process: a phage hangs on the surface of a bacterium cell, injects its 
DNA in it and then forces the bacterium to yield its own replicas inside the 
cell. When the cell is full, it bursts, releasing a huge quantity of copies of the 
infecting phage. The efficiency of the attack (i.e. the probability of success 
of the infection), depends on the couple bacterium-phage. One of the other 
characteristics of bacterial and phages populations are their high variabilities. 
They frequently mutate, creating new populations with new properties. Such 
a system has two dynamics to be taken into account: the dynamics of the 
populations of bacteria and phages and the meta-dynamics of evolution geared 
by mutations and extinctions. It is without any doubt a model which is not 
in the scope of the classical theory of dynamical systems, but our formalism 
applies well to it. 

Dynamical level: it is made up of a set of ordinary differential equations. 
Let us consider the following system 



with S the concentration of nutrient, Bi the different bacteria strains with 
< i < Ub and Pj and Ikj the different phages strains with < j < np. 
To each population Bi and Pj is associated a binary string G {0, 1}" and 
G {0, 1}". These binary strings code for the properties of attack (for the 
phages) or defence (for the bacteria) facing the infection (see [?] for the bi- 
ological discussion). The populations are then characterized by the different 
concentrations and the two lists of bit strings. The set L of all the pairs of 
lists is taken as our second index (the one called J in definition [1]). 

The model is described by a modified version of Mosekilde equations (see 
[?]). This set of ordinary differential equations describes the interactions of 
bacterial populations Bi and phage populations Pj in a chemostat. Bi and Pj 
also symbolize the concentration of these populations and S the concentration 
of the nutrient. The process of infection of bacteria by phages is modelled by 
three infection stages Iij, I2J and I^j. One associates with the ecosystem 
dl]) the state space Xi^nt+iup = IR^"*""''^^"'''. Over a period of time without 
appearance or disappearance of any strain of bacteria or phages, the dynamical 
evolution of the system (jl]) is modelled by the evolution of the concentrations 
of S, the different bacteria Bi and phages Pj, and three infection stages Ikj- 



{S, {-Bi}o< i<nh 1 



{Pj, hj, hj, h,j}o<j<nJ 



(4) 
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With the formahsm of definition [H that means that Tq = M and for a fixed 
/ G L the transition function ipi^nb+inp,i is the integration of the following set 
of differential equations 



= ^ - i?. E au,,P, - pB, 



dl_ 

dt 
dt 
dt 



dPj ^ oPh 
dt T 



i=l 

l{h,3 - hi) - Phi 

' nt rip 3 



dS 
dt 



P, E + E E 

\i=l i=i fc=i ^ 

"6 



wySBj 

K+S 



(5) 



with < z < -n,;,, < J < rip, p the rate of dilution (p = 0.0045 min~^), 
Kand z/ respectively the saturation term and the growth from the Monod 
equation (k = 10pg.ml~^, u = 0.024 min~^), a the theoretical adsorption 
constant depending on phage and bacterium size {a = 10~^ml. min~^), t 
a time constant (r = 30 min), (3 the number of copies of phage j released 
during the burst of the infected bacterial cell (/? = 100), a the continuous 
supply of substrate (a = 10pg.ml~^), 7 the amount of nutrient consumed in 
each cellular division (7 = O.Oln^f) and finally Uij the probability of infection 
of Bi by Pj which depends on the similarity between bit string (attached to 
bacterium population Bi) and (attached to phage population Pj) as follows 

= V — ~) 

with He the size of the binary string. 

Here, we are not interested in the identification of the biological dynamical 
level which can be found in [?]. With the given size of a binary string Uc, 
there exists a finite number of possible dynamical systems, here differential 
equations. Indeed, the number of possible different populations of bacteria is 
equal to the number of parts of the set of binary strings of size Uc- So, there are 
2"^""" possible populations of bacteria. For the same reason, one deduces that 
the possible number of populations of phages is the same and thus the total 
number of possible state spaces is 2"^"" x 2^"'' = 2'^""'^^. So, with the notation of 
definition [H it means that card (J) = 2'^""'^^. If we take for example Uc = 10, 
one has 2^°^^ ~ 3.6 x 10^°^ possible state spaces. Theoretically, one can consider 
a 3.6 X 10'^°^ dimensional space to embed the system ([5]). Nevertheless, when it 
comes to numerical simulations, such a big system is impossible to deal with. 
One thus understands the need of an adaptive system to describe the system 
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Meta-dynamical level (meta-dynamical adaptive system point of 
view): this is the main difference with the model of Mosekilde given in [?, 
Chapter A] which is not evolutive. 

Proposition 11 Consider a small interval of time At, the adaptive changes 
of the system ^ are given by the following mechanisms 

A (t)'' 

p^{t + At, t, 6, (p2+n,+4np,h) " (^1, • • • , 4, (pi+ks+n,+4np,h) = ^^^^^ 

and 

'^{t + At, t, (5, (p2+nt,+4np,h , ^s) = (/?2-E4+nt,+4np,«2 (7) 

where 



(1) formula ^ is the probability that the species S (a bacterium B or a 
phage P) gives birth to k > 1 mutante strains 6i,. . . ,6k on [t,t + At[ 
with A (t) = ^Pe, Pe the probability that a small group of mutant species 
of size 6e gradually replacing the species of the parent population 6. (pij 
is the function defined by equations ^ for j = li,l2 and ks = k if 6 is a 
bacterium and kg = Ak if 6 is a phage, 

(2) formula ^ gives the determinist rule of the extinction of the species S 
which depends on a given threshold rris, es = 1 if 5 is a bacterium and 
es = 4: if 6 is a phage. 

p'^ is the stochastic transition function of definition [7] which governs the meta- 
dynamical rule of the appearance. is the transition function of definition [T] 
which governs the meta-dynamical rule of the extinction. 

Proof. Consider the probability PBiitjk) that the population Bi gives birth 
to k mutante strain on [t, t + At[ and the probability pe that a small group of 
mutante bacteria of size Be is gradually replacing the bacteria of the parent 
population Bi. Such a reasoning gives a binomial probability for 

with n(t) = When n(t) is large, one may approximate the binomial 

probability by the Poisson probability 

p^,(t,A:)^e-^wA^ (8) 

with A (t) = ^^Pe- Suppose that the birth of all the populations to a Ham- 
ming distance of one (only one bit is different) is equiprobable. The birth of a 
mutant strain results in the change of a group of bacteria (which is a part of a 
parent population) of size -Bg, in a population group with new characteristics. 
If this population group already exists, the mutant population is added to 
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it. The same mechanism governs the modelhng of the phage mutation. For 
the extension, there exists a threshold below which the extinction of the 
population is certain. Every population under a given threshold (different for 
bacteria and phages mp) is removed from the system. Thus, the extinction 
is determinist. ■ 

We have defined the macroscopic birth of a mutant population as an event 
occurring on an interval [t,t + At[. There exists a time set Ti where the 
system commutes. This commutation depends on the state of the system 
and on the coefficient Pb pp respectively coefficient of the Poisson law 
of the bacterium and the phage. In order to simplify the model and to make 
it computable, one may suppose that 

Ti = to + i^t 

even if in practice, the mutations have no reasons to be periodically defined. 
This allows to define, with the previous notations, the meta-dynamical rule 

p^{to + iAt, 6, V?2+n6+4np,h) " {^U ■■■,6k, ipi+ks+nt+4np,l2) = ^'^^^^ 

Then, at each step At there are four possible commutations: 

• birth of a new bacterial strain: a variable -Bnj,+i is added, the dimension of 
the system ([5]) increases by one, 

• birth of a new phagical strain: four variables are added: Pnp+i, h,np+i, h,np+i 
and h^np+i, the dimension of the system ([5]) increases by four, 

• extinction of a bacterial strain: the concerned variable is removed, the di- 
mension of the system ([5]) decreases by one, 

• extinction of phagical strain: variables of the concerned phage are removed, 
the dimension of the system ([5]) decreases by four. 

The different possible state spaces resulting from a commutation are given on 
figure m 

If more than one event occurs at each moment t E Ti, one composes the pos- 
sible change of state spaces (for example, the extinction of a phage combined 
with the birth of a bacterial strain decreases the dimension of the system ([5]) 
by 4- 1 = 3). 

Meta-dynamical level (adaptive differential equation point of view) : 

we have described the biological model as a stochastic meta-dynamical adap- 
tive system by using definition [71 This modelling corresponds to a stochastic 
view of the system (j4]) where the space of all stochastic realizations is infi- 
nite. There exists a family of points {tjjjgp^ where the system ([5]) commutes. 
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Appearance of a new 
strain of bacteria 




Fig. 4. Evolution of the state space with Ub and Up respectively the number of 
bacterial and phage strains before the transition. 

We have supposed that this commutation is given by ti = to + i^t. Then, 
at each step At the system may commute. To see a possible evolution of the 
system we choose a stochastic realization at each point t, (only for the ap- 
pearance of a strain because the extinction of a strain is determinist) g{ti,yi) 
where = lim Xi-i (t) with Xi-i(t) the solution of the system ([5]) on [tj_i,ti[. 



It describes the determinist evolution of the system 
a stochastic realization is a function 



By extension, for us 
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X 



ndim{X) 



ndim(X) 



defined at least on {iti,yi)}-^^. For a given realization g, the system (jlj) may 
be modelled by an adaptive differential equation whose equations are given 
by the system ([5]). This modelling belongs to the variable dimensional space 
j^ni,+4np+i _j_ _j_ ]^ ciepends on g and follows the rule given on figure 

m A detailed study of this biological model with implementation can be found 
in m. 



In this example, one sees that the framework of the adaptive differential equa- 
tions corresponds to the case of a transition function defined by a set of dif- 
ferential equations and a meta-dynamical time reduced to a discrete subset of 



6 Conclusion 



We started our research with a biological system whose dynamics changes 
in different dimensions. Seeing that there was no mathematical framework 
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to describe such a system, we have developed a mathematical tool following 
Kalman's dynamical system called meta-dynamical adaptive system which 
was appropriate to give a constructive algorithm for some fractals. It was 
also adapted to describe and analyze our biological system. However, we have 
found our tool too general and we decided to develop a special tool for dif- 
ferential equations. The new system called "adaptive differential equation" 
is not a succession of differential equations because the initial condition of 
each system depends on the previous system and gives the new dimension 
of the following system. This model allows to describe a system of changing 
differential equations, in particular the stochastic realization of a stochastic 
meta-dynamical adaptive system. The last tool we use is the "variable dimen- 
sion space" . This new kind of space links the notion of space and dimension in 
a changing dynamics. We think that our work will contribute to understand 
the huge number of complex systems where the espace of exploration is too 
big to be investigate with classical means. 
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